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SDMMAEY 

Tlie mean and turbulent momentum terms In fully developed turbulent 
pipe flow were experimentally evaluated. The terms of the longitudinal- 
direction momentum eq.uation, obtained from the Eeynolds eq.uations of 
turbulent fluid motion, were experimentally evaluated in a 4 -inch-diameter 
pipe from total- eind static-pressure data and hot-wire anemometer surveys. 
Measurement of the terms appearing in the radial-direction momentum equa- 
tion iudicates the existance of terms as large as or larger than the terms 
of the' longitudinal-direction momentum eq.uaticn. 

Analysis of the turbulent stress tensor shows that the direction of 
principal stress was oriented nearly parallel to the wall in the region 
near the wall. Variation with Eeynolds number of the longitudinal tur- 
bulent intensity' at the center of the pipe indicates that the intensity 
at the center was of a universal nature. 

A direct comparison was made with turbulence measurements obtained 
using the constant -current and constant-temperature systems of hot-wire 
einemometry. The two systems agree well writhin the experimental accuracy 
of the measurements. The constant-temperature measurements also agree 
with the measurements of turbulence presented by laufer for the seime 
Eeynolds number. 


INTEODtJCTION 

Although turbulence in aerodynamics is of greatest tmpojrfcance in 
connection with boundary layers, basic experimental or theoretical studies 
of the boundary layer entail many difficulties. It has proven more prac- 
tical mathematically and experimentally to investigate the simpler forme 
of turbulent shear flow, such as fully developed pipe and channel flows. 
Detailed study of these simpler flows may lead to a better understanding 
of the behavior of turbulent boundary layers. 
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Eesults of measurements of turtulent flow In a pipe and a channel 
have recently "been reported In references 1 and 2. Eesults of further 
measurements made In a fully developed pipe flow are presented herein. 
These measurements were made not caily as an independent verification of 
the results of reference 1, hut also as a comparison between two distinct 
systems of hot-wire anemometiy. 

Measurements were made in a 4-lnch-diameter pipe, which was origin- 
ally built to facilitate the development and calibration of hot-wire 
eq^uipment at the EACA Lewis laboratory. The Eeynolds- number range (based 
on pipe radius) was from approximately 20,000 to 200,000. The measure- 
ments are similar to those reported in reference 1, but Include only the 
mean and statistical quantities necessary to evaluate the terms in the 
Eeynolds equations of motion. 

The two types of hot-wire anemometry equipment were the constant- 
current (or constant-mean -resistance) system and the constant -temperature 
system. In the past, the constant -current system with compensating a-c 
amplifier has been used almost exclusively for the measurement of tur- 
bulence because a preliminary analysis (ref. 3) indicated that the noise 
level of the constant-temperature system would be prohibitively high. 
However, refinements in the system have reduced the noise level of the 
constant -temperature system to a point comparable with that of the 
constant -current system. The actual recording of turbulence data is, by 
comparison, much simpler with the constant-temperature system. 


APPAEATDS APTD PEOCEDDEE 

Test facility . - The test facility consisted of a 4-inch- inside- 
diameter seamless aluminum pipe 19.5 feet in length connected to an air 
supply through a settling chamber and a series of filters, screens, and 
a honeycomb arrangement, as shoim in figure 1. Eocm air, filtered throu^ 
Fiberglas air filters, was supplied by a constant-speed centrifugal 
blower and was metered by a 4-lnch valve. As experimental work proceeded, 
it was necesseiry to add an increasing number of paper air filters in the 
system to prevent dirt particles from affecting the hot-wire measurements. 
The maximum Eeynolds number was 200,000 when measurements were started, 
but as more filtering was added the maximum Eeynolds number was reduced 
to about 150,000 in the region where fully turbulent flow was obtained. 
Data were taken at Eeynolds numbers of 25,000, 50,000, 100,000, and 
150,000. 

Mean velocity instrumentation . - The mean velocity measurements 
were made with a total-pressure probe having a 0.040-inch outside dia- 
meter (with a wall thickness of approximately 0.003 in.) flattened to 
0.010 inch at the tip so that surveys could be made to a point 0.005 
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inch, from the pipe wall. The prohe was accurately positioned at the wall 
hy use of an electrical contact indicator; surveys were made at a position 
14 inches' from the pipe exit. 

Total pressures were recorded continuously by differential pressure 
transducers, which were referenced to atmosuheric pressure. The two 
transducers employed had full-scale pressure ranges of iO.05 pound per 
square inch and i£).15 pound per square inch, with their outputs being 
reco 2 rded on a 2-millivolt full-scale recording potentiometer. Ih the 
systems employed it was possible to detect pressure differences of 
the order of 0.004 pound per square foot and 0.01 pound per square foot, 
respectively, for the two Instruments. A more elaborate description of 
the transducing systems is given in reference 4. 

It has been found experimentally (ref. 5) that total-pressure-probe 
corrections are necessary at very low velocities. In the course of 
measurements in the pipe, it was possible to determine a correction for 
the total-pressure-probe error by comparing measurements of the velocities 
obtained from hot-wire and total-pressure probea at the center of the 
pipe for very low Beynolds numbers. The correction resulted in an in- 
crease in measured velocity of 15 percent at 10 feet per second. There- 
fore, the low velocities in the region near the pipe wall were corrected 
for the probe error. No corrections "were made for probe displacement 
effects due to total-pressure gradients. 

Static-pressure instrumentation . - The pipe •was equipped -with wall 
static taps along the top. The taps were equally spaced 2 inches apart 
in the region where surveys were taJcen and 4 inches apart upstream. The 
pressure distributions for Eeynolds numbers of 25,000 and 50,000 are 
shown in figure 2. A slight bloctoage effect was observed when probes 
were mounted in the pipe, eind corrections were made in e'valuatlng the 
meein velocity and the wall shearing stress as the probe position was 
vaided. These connections were the result of a max lTTimn increase in 
static pressure of 8 to 10 percent as the probe was moved from the wall 
to the center of the pipe. 

Hot-wire turbulence-measuring equipment . - Turbulent fluctuation 
measurements have been made with two Independent systems of hot-wire 
anemometry. Both the conetant-mean-ourrent and constant-instantaneous- 
temperature systems were employed in the measurement of the turbulent 
velocity fluctuations, mainly in an attempt to determine the validity of 
employing the constant-temperature system for boundary layer research. 

The constant- current system (fig. 3) was basically that employed by 
the National Bureau of Standards (ref. 6). The circuit consists of a 
noninductlve-type Wheatstone bridge and an a-c amplifier equipped •with a 
variable resistance-capacitance network to compensate for theimal lag of 
a hot wire. The thermocouple RF meter circuit used to obtain a measure 
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of the mean square of the fluctuating voltage is also included in figure 
3. Figure 4 shows the gain with freq.uency of the complete anemometer 
system with a wire installed in the bridge, both with amplifier com- 
pensation and without compensation. 

The constant-temperature system is that described in reference 7 
with minor changes made to reduce the over-all noise level of the system. 
The output of the constant-temperature anemometer was read on the average- 
sq.uare computers described and evaluated in reference 7. An outline of 
the connections used to obtain the data necessary to evaluate the statis- 
tical averages of the turbulent q.uantities is included in reference 4. 

A wiring diagram of the sum and difference circuit employed in the 

measurements of and is shown in figure 5. This sum and 

difference circuit has proven superior to the circuit reported in ref- 
erence 7, for the noise level and thus the difference Tnlninimn has been 
greatly reduced. The pfirticular instrument shown in figure 5 was de- 
veloped at the Lewis laboratory, but the general idea is to be found in 
several commercial computing instruments. 

The hot-wire probes are the same as described in reference 4. As 
discussed therein, two possible methods of wire calibration were avail- 
able. Initially, the calibration procedure was to place the wire at the 
pipe center line and record the heat loss from the wire as the pipe mass 
flow was varied - the mass flow PU being determined by Independent 
measurements with a total- pressure probe. The second method consisted 
of a continuous calibration of the wire as it was operating in the pipe. 
Talues of -\/OT at each r disteince obtained from mean velocity profile 
measurements with a total-piassure probe plotted against the wire heat 
loss recorded when the fluctuation measurements were taken results in a 
continuous calibration as a turbulent profile is taken. The continuous 
(second) method of calibration was used, since dust particles in the 
air stream caused discontinuous changes in the calibration curve. 

The method of least squeires was used to determine the elopes of the 
calibration curves. Because of the aforementioned discontinuous changes 
in wire calibiation, it was often necessary to decide arbitrarily which 
particular set of points belonged to a given portion of the calibration 
curve. Cross checks between measurements (measurements at the pipe 

center line with Reynolds number and measurements evaluated from 

both the single and x-wire p37obes) were therefore employed in determining 
the final calibration. It is believed that the maximum error due to wire 
calibration was no greater than ±10 percent. A typical continuous-type 
calibration curve is shown in figure 6. 

Some deviation from the linear relation between heat loss and the 
sq.uare root of mass flow was noted for the high velocities (180 ft/sec 
and greater) . 
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DECDBSION OJ EESDLTS 


Completely descsri'bliig tlie turbulent flow field necessitates knowing 
not only the time averages of the turbulent quantities appearing in the 
momentum and energy equations, but also all hi^er-order terms that mi^t 
be obtained by taking hi^er moments of the momentum equation. At pre- 
sent it is neither practical nor experimentally possible to obtain all 
such terms; thus esperlmental work has been directed mainly towaord mea- 
surements of only those terms appearing in the momentum equation' and its 
first moment. The present measurements are limited almost entirely to 
the momentum teims. All terms of the equations of motion, simplified 
for fully developed pipe flow, have been measured or can "be evaluated 
for Reynolds numbers of 25,000 and 50,000. Supplemental data at Reynolds 
numbers of 100,000 and 150,000 are also presented. 

Mean velocity distributions . - Mean velocity profiles for the four 
Reynolds numbers where turbulence measurements were taken (Re = 25,000, 
50,000, 100,000, and 150,000) are presented in figure 7(a). An enlarged 
portion of the velocity profiles near the wall is shown in figure 7(b). 


The elope of the velocity profiles at the wall as predicted from 
dp/dx is also Included in figure 7(b). Jigure 8 shows a plot of the 

viscous shearing stress coefficient ^V/^7 near the wall. The eiis- 

ience of a nonvanishing second derivative of the velocity distribution 
at the wall (as indicated in fig. 8) was predicted by the analysis of Pai 
(ref. 8), which indicates that a nonlinear velocity distribution is con- 
sistent with the equations of motion and the boundary conditians. 


Evaluation of terns in z-dlrectlon momentum equation . - The x- 
dlrection equation of motion for fully developed pipe flow may be reduced 
to (ref. 1) 


uv 





( 1 ) 


For a major portion of the flow field liv will be the Important shear 
term present, since only near the wall (fig. 8) is the viscous shear 
stress large. The distributions of Puv, in coefficient form, across 
the pipe are shown in figure 9. The right side of equation (1) is plotted 
nondlmensionally as the dashed curves of figure 9 and may be compared 
directly with the measured values of puv. Althou^ there is consider- 
able scatter in the uv measurements, it is evident that the data satis- 
fy equation (l) for the range of Reynolds number investigated. 


For the two higher Reynolds numbers (figs. 9(c) and (d)), it was 
necessary to correct the output of one of the hot wires (approximately 
40 percent) so that a zero shear stress would be indicated at the center 
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of the pipe. (The error apparently was due to dsTlation of the eingle of 
one wire^ as hoth profiles were taien with the same prohe and the same 
trouble was encountered.) Uo corrections were necessary for the prohes 
used at the two lower Reynolds numbers. The points very close to the 
wall are also included, but it is evident that the presence of the wall 
affected the heat loss from the wires. 

The measurements at a Reynolds number of 25,000 are compared with 
those reported in reference 1 for a 10-inch-diameter pipe, which were 
obtained by a method slightly different from that used in this report (see 
ref . 4 for method employed herein) in that the correlation coefficient 

was measured directly, and uv was calculated from the 

/iff 

coefficient and independent measurements of 'V u and ty v^ . The method 
of reference 1 appears to have resulted in less random scatter in the 
values of Puv, but it does require a wire matching accomplished using 
the assumption that uv = 0 at the center of the flow. 

Evaluation of terms in r-direction momentum equation . - The r- 
direction equation of motion has in the past received little or no atten- 
tion in either esjiertmental or theoretical Investigations, since for 
laminar shear flows the terms of this equation were found to be negligi- 
ble compared with those of the s-direction eq.uation. However, for tur- 
bulent shear flows, it is found that the turbulent contribution to the 
r-direction momentum results in terms as large as or larger than those 
appearing in the x-direction equation. It may be sho^jn rigorously, how- 
ever, that for fully developed turbulent pipe flow the i- and r- 
directlon equations of motion are Independent (ref. 1). 

The r-direction equation of motion can be reduced to the form 
(ref. 6) 




- V^ ^ Pw " P 
dr = — - — r 


( 2 ) 


Evaluation of all terms in equation (2) requires the measurement of v^ 

o 

6uid w along with the variation of the static pressure across the pipe. 

The measurement of v^ and w^ is readily accomplished by hot wires. 

The static-pressure variation across the pipe is no greater th£m a few 
thousandths of an inch of water; thus no attempt was made to measure it. 


The variations of and across the pipe for Re = 

25,000 and Re = 50,000 are shown in figures 10(a) and (b) and 11(a) and 
(b), respectively. It can be seen that the two velocity components are 

of the same order of magnitude, with the values of rising to 
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sli^tly higlier values as the -(rail is approached. The va.lues obtained 
from the faired curves of figures 10 and 11 have been used to evaluate 
the terms on the left side of equation (2), and the resulting nondimcjn- 
sionalized distributions are presented in figures 12(a) and (b) . Talues 
of the right member of equation (2), obtained by subtraction of the two 
terms in the left member ^ are also plotted and show that, although the 
absolute magnitude of the pressure change across the pipe is small, 
the timisverse pressure gradient is quite large, particularly near the 
trail. As a comparison of the orders of magnitude of the terms of the 
X and r momentum equations, the ri^t side of equation (1), as pre- 
sented in figure 9, has also been plotted in figure 12. The existence 
of large terms in the y-direction equation of motion for boundary layers 
has also been discudsed in reference 4. 


Turbulent stress tensor . - A study of the momentum distribution 
cannot be limited wholly to the evaluation of the time-average equaticnc 
of motion, since in these equations the averaging process has eliminated 
some of the momentum terms present in the flow field. Eeynolds' (ref. 9) 
Investigation of turbulent flow resulted in the identification of certain 
virtual stress terms eirising from the turbulence, which can be shown to 
form a turbulent stress tensor 





u^ uv uw 
vu v^ vw 
wu wv w^ 


( 3 ) 


For the fully developed turbulent pipe flow it is possible to show that 

wv = vw=wu=uw = 0 and uv = vu 


so that only four of the terms are independent. 

The only tom appearing in the stress tensor which does not also 

O 

appea r in scane form in the equations of motion is Pu . The measurements 

of for the four Eeynolds numbers investigated are presented in 

figure 13. It may be seen that the values of u^ are by far the largest 
terns in the stress tensor (eq. (3)), with magnitudes as great as four 
times that of any of the other terms for the region close to the wall. 

Near the center the value of u^ approaches that of v^, so that the 
conditions of isotropy are nearly satisfied at the center of the pipe. 

A second factor of Importance in regard to the turbulent stress 
tensor is the direction along which the maximum stress acts. For the 
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■tui'bulen't strees ’tensor, tlie angle p 'be’tveen the principal axis of the 
tensor and the direction of ’mean flow is given hy the equation (ref. 10) 


tan.2p 


2ut 



( 4 ) 


The ■values calculated for P are shown In figure 14. The evaluation of 

P near the center of the pipe is not accurate because equation is 

2 2 

extremely sensitive to the measured values in this region since u v 
and uv -> 0 at r = 0. As the -wall is approached, the angle P decreases 
so that the direction of principal stress is alined nearly parallel to 
the flow and the -wall. 


Beynolds number effects. - The effect of Eeynolds number on turbu- 
lent intensity has been deteimined from experimental measurements of the 
turbulent intensities at the pipe center. Figures 15 and 16 show the 

■variation of the -two 'turbulent velocity components and 

along the pipe center -with Eeynolds number. The points are also compared 
with da'ta from references 1, 2, and 11. Ee suits in figure 15 indicate that 
the -variation of the longitudinal component -with Ee can be considered 
universal within the limits of experimental accuracy. It is noted in 
reference 12 that for a fully developed channel flow the Intensity at the 
center -was Independent of any upstream turbulence, thus showing that the 
turbulence has reached a statistical eq.uilibrium state which will not 
change unless the symmetry of the flow is disturbed. The following em- 
pirical eq.uation is found to fit the data of figure 15: 


(is: 

I u_ i 


= 0.144 Ee 


-0.146 


° /at r=0 


(5) 


A summEiry of the variation of the longitudinal turbulent intensity 

distribution is sho-wn in figure 17. Figure 18 shows the profile 

presented in reference 1 for Ee = 250,000 compared with a distribution 
for Ee = 200,000 measured in the 4-inch pipe. 

Comparison of cons’bant- current and constant-temperature hot-’wire 
anemometer systems . - A direct comparison of measurements made using 
the two types of hot-’wire anemometer systems employed is shown in figure 
15. The feasibility of using the constant-temperature system for measure- 
ment of turbulence intensities encountered in -burbulent shear flo-ws has 
been q.uestioned (ref. 3), since the noise level of the system was con- 
sidered too high. However, for the turbulence levels encountered in the 
pipe, cons'cant-temperature measurements agree -with constant -current 
measurements ■within the range of experimen'tal accuracy. All distribution 
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measurements across the pipe reported herein vere taken with the constant- 
temperature instruments, since only a fraction of the time per measure- 
jnent was req.uired as compared with that needed in constant-current 
measurements. Since the measurements reported hy Laufer (ref. 1) were 
perfoimed with a constant-current system, the comparison of profiles 
(figs. 9(a), 10(a), 11(a), and 13(a)) also offers a check of the agree- 
ment "between the two systems. 


Comparison of pipe and channel flow . - A brief comparison "between 
the flow in a pipe and that in a channel was made. Measurements of the 



distribution in a channel for Be = 61,000 reported in reference 2 


are shown in figure 19. The distribution of ^u^ measured for the same 
Eeynoldd number in the 4-lnch pipe is also Included in figure 19. It 
may be seen bhat the distributions for the pipe and the channel are 
similar. 


Specti»um measurements . - The spectra of the longitudinal turbulent 
Intensity for Be = 25,000 at several r-distances obtained by use of a 
harmonic wave analyzer are shown in figure 20. The function J(n) was 
defined as 


F(n) 


^(^n) 


u^ dn 


( 9 ) 


where dl 
dn 


.^u^^ is the contributlan to u^ from freq.uencies n - ^ and 

n + — , and dn is the effective band width of "che measuring analyzer 
2 

(4.4 cps) . Figure 20 shows a redistribution in the hi^-freq.uency range 
for the spectrum distribution very near the wall. 


CONCIODING EEMAEKB 

Analysis of the measured momentum terms in fully developed turbulent 
pipe flow indicates that the balance of the longitudinal-direction momen- 
tum terms as predicted from the Eeynolds equation of motion can be ex- 
perimentally verified. Evaluaticm of the terms in the radial-direction 
momentum balance shows the existence of a pressure gradient across the 
pipe. The terms of the r-direction momentum equation for turbulent pipe 
flow are comparable in magnitude with those of the x-directlon equation; 
however, the two equations are independent. 

Evaluation of the turbulent stress tensor shows the mean square 
longitudinal -velocity to be by far the largest term near the wall, with 
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tlie other three terms heing of roughly eq.ua! magnitude. The calculation 
of the angle between the principal axis of the stress tensor and the x- 
direction shows the axis to he oriented nearly parallel to the wall in 
the region near the wall. 


A study of the Eeynolds number effects on the longitudinal turbulent 
Intensity at the center of the pip® indicates that the intensities are 

of a universal nature. Comparison of and at the center 

indicates the flow to be nearly isotropic. 


The longitudinal turbulent intensity at the center of the pipe could 
be represented by the empirical relation 



0.144 Be 


-0.146 


where 'ip’ is the root mean sq.uare of the longitudinal turbulent 
velocity, TJq is the mean velocity at the center of the pipe, and Be 
is the Beynolds number. 


Direct comparison of the constant-current and constant -temperature 
systeias of hot-wire measurements of turbulent q.uantlties showed agree- 
ment within the limits of experimental accuracy. The constant -temperature 
measurements at Be = 25,000 are found to agree within the limits of ex- 
perimental accuracy with the measurements presented by Laufer for the 
same Beynolds number. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 27, 1954 
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APPERDIX - SIMBOLS 

a radius of pipe 

h. channel half -width 

i mean current of hot wire 

p mean static pressure 

p^^ mean static pressure at wall 

E operating resistance of hot wire 

E„ resistance of hot wire at reference condition 

Ct 

Be Eeynolds number based on channel half -width or pipe radius and- 

velocity at center of flow, pU^a/n or pUgh/ja 

r radial distance from pipe center 

r' (a - r) 

Tm turbulent stress tensor 

U mean longitudinal velocity 

TJq mean longitudinal velocity at center of pipe or channel 

TX^ shear stress velocity, -y%^/p 

u,v,w longitudinal, radial, and tangential components of turbulent 
velocity 

uv Eeynolds turbulent shesir stress in zr-plane 

uw Eeynolds turbulent shear stress in x-tangential plane 

vw Eeynolds turbulent shear stress in r-tangential plane 

X direction of mean flow 

P angle betireen principal eixis of turbulent stress tensor and 

direction of mean flow 

p viscosity of air 

V kinematic viscosity of air 



12 


HACA TN 3266 


p density of air 

shearing stress at wall 
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Figure 6. - Hot-wire calibration curve. Constant-temperature operation; 
slope calculated "by method of least sq.uEires neglecting points near wall. 
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Fi^cure 9. - CoTrtdniiftd.. pla*tri*birfcioTi of tur'buloirb sb.e€Lr fl*bresa coefficient across pipe. 
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Cb) Reynolds nunber, 50,000. 

Figure 10. - Radial turbulent Intensity distribution across pipe. 
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(a) Reynolds ntnaber, 25,000. 



(b) Reynolds nunber, 50,000. 

Figure 11. - Tangential torbulent intensity distribution across pipe. 










(a) Reynolds number, 25,000. 

Figure 13. - Longitudinal turtulent Intensity distribution across pipe. 
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?lgure 15. 

- Longitudinal turhulont Intensity at centor of plpa. 
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Figure 18. - Longitudinal turtiulent intensity distriLutlon across pipe for large Reynolds numbers 
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Figure EO. - Spectra of longitudinal turlulent energy. Reynolds 
number, 25,000. 


NACA - Langley Field, V*. 



